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Abstract — ^Multistatic  radars  utilize  multiple  transmitter  and 
receiver  sites  to  provide  several  different  monostatic  and  bistatic 
channels  of  observation.  Multistatic  passive  and  active  systems 
can  offer  many  advantages  in  terms  of  coverage  and  accuracy  in 
the  estimation  of  target  signal  parameters  but  unfortunately 
their  performances  are  heavily  sensitive  to  the  position  of 
receivers  (RX)  and  transmitters  (TX)  with  respect  to  the  target 
trajectory.  As  known,  geometry  factors  play  an  important  role 
in  the  shape  of  the  ambiguity  function  (AF)  which  is  often  used 
to  measure  the  possible  global  resolution  and  large  error 
properties  of  the  target  parameters  estimates.  Exploiting  the 
relation  between  the  ambiguity  function  and  the  Cramer-Rao 
lower  bound  (CRLB),  in  this  work  we  propose  an  algorithm  for 
choosing  in  a  multistatic  scenario,  along  the  trajectory  of  the 
tracked  target,  the  pair  TX-RX  with  the  best  asymptotic 
performance  calculated  in  terms  of  CRLB  on  estimation 
accuracy. 

L  Introduction 

Recently,  great  interest  has  been  devoted  to  systems 
making  use  of  illuminators  of  opportunity,  such  as  broadcast 
or  communications  signals,  tracking  targets  by  range  and 
Doppler  information.  These  techniques  are  known  as  Passive 
Coherent  Location  (PCL),  and  have  the  advantage  that  the 
receivers  do  not  need  any  transmitter  hardware  of  their  own, 
and  are  completely  passive,  and  hence  undetectable.  Among 
all  the  transmitters  of  opportunity  available  in  the 
environment,  broadcast  transmitters  represent  some  of  the 
most  attractive  for  surveillance  purposes,  owing  to  their  high 
powers  and  excellent  coverage. 

Furthermore,  Passive  Coherent  Location  techniques  can  be 
extended  to  Multistatic  radar  systems  where  multiple 
transmitter  and  receiver  sites  are  used  to  provide  several 
different  bistatic  channels  of  observation,  leading  to  an 
increase  in  the  information  on  a  particular  area  of  surveillance. 
The  information  gain  obtained  through  this  spatial  diversity 
can  give  rise  to  a  number  of  advantages  in  typical  radar 
functions,  such  as  detection,  parameter  estimation,  tracking 
and  identification.  The  performance  of  each  bistatic  channel 
heavily  depends  upon  the  geometry  of  the  scenario  and  the 
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position  of  the  target  with  respect  to  each  receiver  and 
transmitter. 

In  this  paper  we  approach  the  problem  of  optimally 
selecting  the  transmitter-receiver  (TX-RX)  pair,  through 
calculation  of  the  ambiguity  function  (AF)  and  of  the 
Cramer-Rao  lower  bounds  (CRLBs)  for  the  bistatic  geometry 
of  each  TX-RX  pair.  The  best  pair  is  defined  as  that  exhibiting 
the  lowest  bistatic  CRLB  for  the  target  velocity  or  range. 
These  results  can  be  used  for  the  dynamical  selection  of  the 
TX-RX  signals  for  the  tracking  of  a  radar  target  moving  along 
a  trajectory  in  a  multistatic  scenario.  The  multistatic  scenario 
analyzed  in  this  work  is  composed  by  a  transmitter  and  four 
receivers.  The  results  reported  here  have  been  obtained 
considering  a  single  frequency  modulated  (FM)  commercial 
radio  station  as  transmitter  of  opportunity  [1]. 

II.  Bistatic  Geometry 

Before  starting  our  study,  it  is  necessary  to  describe  the 
coordinate  system  used  to  represent  a  bistatic  radar  geometry. 
Figure  1  shows  the  coordinate  system  and  its  parameters.  The 
positions  of  the  TX,  RX  and  target  are  generic.  Considering  an 
ordinary  Cartesian  grid,  the  TX  is  located  at  point  7x,  whose 
coordinates  are  (x^,  y^),  the  RX  is  at  point  Rx  in  {xr,  yR)  and 
the  target  is  located  at  point  B,  whose  coordinates  are  (x,  y). 

The  triangle  formed  by  the  transmitter,  the  receiver  and  the 
target  is  called  the  bistatic  triangle.  As  shown  in  Figure  1,  the 
sides  of  the  bistatic  triangle  are  Rr,  Rr  and  L,  where  Rr  is  the 
range  from  transmitter  to  target,  Rr  is  the  range  from  receiver 
to  target  and  L  is  the  baseline  between  the  transmitter  and  the 
receiver.  The  internal  angles  of  the  bistatic  triangle  that, 
without  lack  of  generality,  are  assumed  to  be  positive,  are  a,  S 
and  y.  In  particular,  the  bistatic  angle  S  is  the  angle  at  the  apex 
of  the  bistatic  triangle,  at  the  vertex  which  represents  the 
target.  Assuming  that  the  coordinates  of  the  transmitter,  the 
receiver  and  the  target  are  known,  it  is  possible  to  calculate  all 
the  parameters  of  the  bistatic  triangle.  Or  and  Or  are  the  look 
angle  of  the  transmitter  and  the  look  angle  of  the  receiver, 
respectively,  they  are  measured  positive  clockwise  from  the 
vector  normal  to  the  baseline  pointing  forward  the  target. 
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Figure  1 .  Bistatic  geometry 


From  Figure  1,  we  have  that:  Q^TiH-a,  6R=y-7il2, 
S=7i-a-y=0T-0R,  and  from  the  cosine  law  we  obtain: 
RT=RR+L^-2RRLsmOR,  which  gives  the  range  from 
transmitter  to  target  Rr,  as  a  function  of  the  range  from 
receiver  to  target  Rr  and  the  look  angle  of  the  receiver  Or.  In 

Figure  1  we  plot  also  the  target  velocity  vector  V  ;  ^  is  the 
angle  between  the  target  velocity  vector  and  the  bistatic 
bisector,  which  is  measured  in  a  positive  clockwise  direction 
from  the  bisector.  In  particular  the  bistatic  bisector  is 
represented  by  the  vector  BI ,  where  /  is  the  incenter  of  the 
bistatic  triangle,  whose  coordinates  are  (x/,y/). 

In  the  bistatic  geometries,  an  important  parameter  is  the 
radial  velocity  F^,  which  is  the  target  velocity  component 
along  the  bistatic  bisector.  From  the  observation  of  Figure  1, 

we  obtain  =  V  ■  =  |f|  cos  (j) . 


III.  Bistatic  Ambiguity  Function 

In  our  study  we  suppose  that  the  complex  envelope  of  the 
signal  transmitted  by  the  FM  commercial  radio  station  is  the 
unitary  power  pulse  given  by: 


u{t)  =  < 


1 

VF 


^jp?,m{27tfQt+(p) 


0 


elsewhere 


(1) 


that  is  a  signal  which  instantaneous  frequency  is  a  sinusoidal 
oscillation  e  [2-3].  In  particular,  T  is  the  observation  time,  is 
the  modulation  index  and  \/fo  is  the  period  of  the 
instantaneous  frequency.  In  other  words,  we  assumed  that, 
during  the  observation  time,  the  modulating  signal 
transmitted  by  a  radio  station  can  be  approximated  by  a 
sinusoidal  oscillation.  This  can  be  justified  considering  that  in 
a  typical  FM  radio,  the  program  content  is  speech  and/or 
music,  which  are  often  modelled  as  periodic  vibrations. 
Moreover,  the  chosen  signal  is  a  mathematically  tractable 


model  that  makes  it  feasible  to  study  the  analyzed  scenario 
rigorously. 

As  widely  known,  the  complex  ambiguity  function  (CAF) 
of  a  pulse  u{t)  represents  the  response  of  a  filter  matched  to  a 
given  finite  energy  signal  when  the  signal  is  received  with  a 
delay  and  a  Doppler  shift  relative  to  the  nominal  values 
Th  and  Vh  expected  by  the  receiver.  Therefore,  the  CAF 
definition  is  [4-5]: 


In  eq.  (2),  and  are  the  actual  delay  and  Doppler 
frequency  of  the  radar  target,  respectively,  and  Th  and  Vh  are 
the  hypothesized  delay  and  frequency.  In  the  monostatic  case 
there  is  a  linear  relation  between  and  and  the  target  range 
position  Ra  and  radial  velocity  Va  ,  i.e.  =2RJc  and 
=  2  V/cIc,  where  fc  is  the  carrier  frequency  of  the  reference 
signal  u{t).  Similar  relations  hold  for  Th  and  V//.  Clearly, 
is  maximum  for  th=  Ta  and  Vh=  v«. 

Based  upon  definition  (2)  we  can  calculate  the  monostatic 
CAF  for  the  signal  u(t)  in  eq.  (1).  In  particular,  after  some 
manipulation,  it  is  possible  to  write 


n  k 

where 


sin(;rF(r-|T|)) 

ttvT 


(4) 


for  |r|  <  T  and  0  elsewhere,  where  we  set  r  =r//-  r^,  v  =V//- 
and  JniP)  is  the  order  Bessel  function  of  the  first  kind. 

As  an  example.  Figures  2-4  show  the  Ambiguity  Function 
(AF),  that  is  the  absolute  value  of  the  CAF,  of  the  signal  in 
eq.  (1)  when /^=15kHz,  P=5,  T=2lfo  and  (p='Kl2.  In  particular, 
the  0-Doppler  cut  is  characterized  by  the  presence  of  two 
secondary  lobes.  Considering  T=klfo  and  keN  ,  the  number 
of  secondary  lobes  is  2(A:-1). 


|X(t.v)| 


1  'T-'"" 


0.8^. 
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Figure  2.  AF  of  u(t),  3D-plot;/o=15kHz,  ^=5,  T=llfo,  (p=nl2. 
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To  refer  eq.  (3)-(4)  to  the  bistatic  geometry  of  Fig.  1  and  to 
obtain  the  expression  of  the  bistatic  ambiguity  function,  we 
must  replace  Th  and  Vh  in  (3)-(4)  with  the  following  formulas 
[6], [7]: 


(5) 


and 


Jc 


=  2—V COS0  |— + 


Rr+L  sin 


2  2^Rl+L^  +2Rj^Lsme^ 


■  (6) 


Unlike  the  ambiguity  function,  which  provides 
information  on  the  global  resolution,  the  CRLBs  are  a  local 
measure  of  estimation  accuracy.  Anyway,  both  can  be  used  to 
asses  the  error  properties  of  the  estimates  of  the  signal 
parameters.  In  [8]  the  author  derived  a  relationship  between 
CRLB  and  AF,  which  has  been  successfully  used  in  the 
analysis  of  passive  and  active  arrays  [5].  In  the  monostatic 
configuration,  [8]  claims  that  for  the  Fisher  Information 
Matrix  (FIM)  the  following  relationship  holds: 


ar" 

d^^{T,v)\ 

dvdr 


dTdv 

dv^ 


T=0,V=0 


(7) 


where  SNR  is  the  signal-to-noise  power  ratio  at  the  receiver. 
Using  the  results  obtained  in  [8],  it  is  also  possible  to 
demonstrate  that 


Figure  4.  AF  of  u(t),  0-Doppler  cut; ^^=1 5kHz,  j3=5,  T=2lfo,  (p=7il2. 
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Figure  5.  AF  of  u(t),  0-delay  eut; /o=15kHz,  j3=5,  T=2lfo,  (p=nl2. 


d^\x{t,v)\ 


dv^ 

d"\X{T,v)\ 

drdv 


-fOO 

=  j  \u{t)f  dt 


(8) 

(9) 


T=0,v=0 


t  u  (t)dt  !>  (10) 


dt 


After  some  algebra,  it  is  possible  to  verify  that: 

=25M?^^^[2;r/„r  +  sin(2;r/„r)cos(2^)]  (11) 


[J^I^=2SNR^ 


(12) 


[Jm],.,  =  -2SNR  [;r/„rcos(;r/„r)  -sin(;r/„r)]  (13) 

No 


The  presence  of  secondary  lobes  is  due  to  the  periodicity  of  ^  x  i  ^ 

the  analyzed  pulse.  Moreover,  for  k>l,  the  AF  has  a  T-  From  eq.  (7)  the  CRLBs  follow: 

periodic  behavior  along  the  delay  axis.  Obviously,  if  Tfo^  1,  CRLB(T^)  =  [j)^(r^,u^)J^^  and  CRLB(u^)  =  [j^(t^,u^)]^^  . 
u(t)  tends  to  a  rectangular  pulse  and  its  CAF  tends  to  Xi(r,v) 
in  eq.  (4). 
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As  an  example,  assuming  T=klfo  (keN)  and  considering 
that  7r^k^/3»sm^(g)),  the  Root-CRLBs  are  approximated  by 


RCRLB(rJ  =  ^CRLBCrJ  = 
RCRLB(i/J  =  ^CRLB(vJ 


1  1 
^2SNR  TtBc 

1  s 

y/lSNR  T 


(14) 

(15) 


where  Bc—ipfo  is  Carson’s  Bandwidth'  of  u{t).  It  is 
interesting  to  observe  that  RCRLB(v^)  is  inversely 
proportional  to  the  time  duration  of  the  reference  signal, 
while  RCRLB(ra)  is  inversely  proportional  to  its  bandwidth. 
From  this  result  we  can  observe  that  the  best  performance  is 
obtained  with  modulating  signals  with  high  spectral  content, 
such  as  rock  music,  and  poorest  performance  is  obtained  with 
slow  varying  modulating  signals,  such  as  speech  signals. 

In  the  bistatic  configuration  we  should  express  the 
ambiguity  function  in  terms  of  the  bistatic  t{Rr,6rX)  and 
v{Rr,Vb,0rJ^),  and  derive  the  AF  with  respect  to  the  useful 
parameters  Rr  and  Vr,  as  follows: 


^"\x{R„Vs)\ 

dR/ 

dF,dR, 


dR,dV, 

^"\x{R,,Vs)\ 

dF/ 


(16) 


RR=RayB=Va 


Using  the  derivative  chain  rule  and  letting  RR=Ra  and  VB=Va, 
after  some  algebra  it  is  possible  to  verify  that  [5]: 


[Js]2,2  “  1 


dr) 
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T  dz  dv  T 
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(17) 

(18) 
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Figure  6.  RCRLBs  as  a  function  of  the  look  angle  Z=50kni,  R7j=30Km, 
Ffi=250m/sec,/o=15kHz,y5=5,  T=2Qlfo,  (^=7r/2,/c=100MHz,  SNR=20dB. 


0  [rad] 


Figure  7.  RCRLBs  as  a  function  of  the  look  angle  Or,  L=50kni,  RR=30Km, 
VB=250m/sQcJo=l5kllz,j3=5,  T=20/fo,  (?J=7r/2,/c^l00MHz,  yA7?=20dB. 


For  L=0  the  bistatic  CRLBs  coincide  with  the  monostatic 
CRLBs.  To  highlight  the  differences  between  monostatic  and 
bistatic  domain.  Figures  6  and  7  show  the  RCRLBs  as  a 
function  of  the  baseline  length  L  and  the  angle  9r. 

It  is  evident  that,  for  all  the  parameter  values  we  tested, 
the  bistatic  RCRLBs  are  higher  than  the  monostatic  RCRLBs. 
The  bistatic  RCRLBs  get  even  worse  when  the  target 
approaches  the  baseline,  that  is  when  Rr<L  and  9r 
approaches  -tt/I.  In  this  case,  the  resulting  delay  is  L/c  and 
the  radial  velocity  is  zero,  therefore  resolution  is  totally  lost 
and  the  RCRLBs  tend  to  infinity. 


Carson's  bandwidth  rule  defines  the  approximate  bandwidth  of  an  angle- 
modulated  signal  and  is  expressed  by  the  relation  Bc=2(Af+B)  where  A/is  the 
peak  frequency  deviation,  and  B  is  the  bandwidth  of  the  modulating  signal. 


IV.  Optimal  Selection  of  the  TX-RX  Pair 

In  our  scenario  we  considered  a  surveillance  map  of 
dimension  100km  and  100km  where  we  placed 
1  transmitter  and  4  receivers. 

In  particular,  we  placed  the  transmitter  at  coordinates 
T={LJ2,  Ly/l)  and  the  receivers  at  coordinates 
LyU),  rF'={LJ1,  1Lyl\\  Rp’'^=(LJl,  V4)  and  R‘''^=(34/4, 
Ly/2).  Therefore,  there  are  4  pairs  of  TX-RX  that  we 
considered  as  independent  bistatic  systems.  The  results  have 
been  obtained  considering  a  non  co-operative  frequency 
modulated  (FM)  commercial  radio  station  as  transmitter  of 
opportunity  and  modelling  the  reference  signal  using  the 
pulse  in  eq.  (1).  The  carrier  frequency  of  the  system  was  fixed 
to  fc=  lOOMHz,  which  falls  within  the  VHF  part  of  the  radio 
spectrum  in  which  FM  is  used  for  broadcasting. 
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In  order  to  fix  the  values  of  fo  and  we  considered  that 
the  signal  emitted  by  a  FM  commercial  radio  station  is  an 
audio  signal  (speech  and/or  music),  therefore  the  frequency  of 
the  modulating  signal  fo  has  to  belong  to  the  range  of  audio 
frequencies,  whose  accepted  standard  range  is  from  20Hz  to 
20kHz.  Moreover,  FM  commercial  radio  stations  use 
bandwidth  of  about  150kHz,  therefore,  according  to  Carson’s 
rule,  we  fixed /o=  15kHz  and  P=5. 

In  the  analyzed  scenario  we  also  considered  an  omni¬ 
directional  antenna  for  the  transmitter  and  directive  antennas 
for  the  receivers.  In  particular,  we  considered  that  the 
antennas  of  the  receivers  have  the  same  power  gain: 


Gm  =  \ 


In  2|6>^+;r/2| 
1k-£  £ 

Ik 

1k-£ 


|^.+V2|<| 

elsewhere 


(20) 


that  is,  the  receiver  antenna  power  gain  is  constant  for  all 
the  directions,  has  a  notch  towards  the  transmitter  {0r=-kI2) 
and  linearly  decreases  when  the  receiver  look  angle  belongs 
to  the  range  [-;r/2-e/2;-;r/2+e/2].  In  our  simulation  we  fixed 
s=kI6.  Wq  choose  this  receiver  antenna  power  gain  because  in 
a  real  case  scenario  the  receiver  antennas  are  designed  in 
order  to  maximize  the  power  ratio  between  the  target  echo 
and  the  Line  of  Sight  co-channel  interference. 

For  each  point  of  the  analyzed  area  and  for  each  of  the  4 
bistatic  channels  we  calculated  RCRLBs  on  the  target  range 
and  velocity  estimation  accuracy.  In  particular,  we  assumed 
that,  in  each  point  of  the  analyzed  area,  the  target  has  a 
velocity  vector  aligned  to  the  x  axis  with  intensity  250  m/sec. 
The  RCRLBs  of  target  range  and  velocity  are  a  function  of 
the  range  from  the  receiver  to  the  target  Rr,  the  baseline  L, 
the  look  angle  of  the  receiver  Or,  the  radial  velocity  Va,  and 
the  Signal-to-Noise  Power  Ratio  (SNR).  All  these  parameters 
depend  on  the  configuration  of  the  bistatic  triangle,  that  is,  on 
the  coordinates  of  target,  TX  and  RX.  Bistatic  geometry  also 
affects  the  received  echo  power,  because  the  path  loss  factor 
is  (RrRr)^.  In  particular,  the  SNR  can  be  written  as 
SNR=G(dR)SNRc(L^^+Lyy/(RTRRf,  where  SNRc  is  a  constant 
parameter.  We  assumed  that  SNRc  =20dB.  That  is,  we 
assumed  that  if  both  the  transmitter  and  the  receiver  are 
located  in  (0,0)  and  the  target  is  located  in  (L^fy),  then 
5’/W?=G(i9^)|dB+20dB. 

Figures  8  and  1 1  are  colour  coded  maps  which  represent 
the  RCRLBs  of  the  target  range  and  the  target  velocity  in 
each  point  of  the  analyzed  area  for  the  4^^  bistatic  system.  In 
particular.  Figure  8  represents  the  RCRLB  of  the  target  range, 
in  dBm  while  Figure  1 1  represents  the  RCRLB  of  the  target 
velocity  in  dBm/sec. 

As  apparent  from  the  results,  the  RCRLBs  of  the  bistatic 
system  are  strongly  related  to  the  bistatic  geometry.  It  is 
evident  that  the  effects  of  geometry  factors  are  more 
prominent  as  the  target  is  in  the  same  direction  of  the 
transmitter,  in  this  case  the  resolution  is  totally  lost.  This  is 
clearly  due  to  the  geometrical  factors  discusses  in  the 
previous  section  and  to  the  choice  of  the  receiver  antenna 
power  gain  in  eq.  (20).  However,  the  effects  of  the  bistatic 


geometry  are  less  prominent  when  the  distance  to  the  target 
increases;  in  this  case  the  bistatic  system  behaves  more  and 
more  as  a  monostatic  system. 


Figure  8.  RCRLB  of  the  target  Range  [dBm],  4^  pair. 


Km 


Figure  9.  Minimum  RCRLB  of  the  target  range  [dBm]. 
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Figure  10.  Optimum  pair  map  for  target  range  estimation. 
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Figure  1 1 .  RCRLB  of  the  target  veloeity  [dBm/see],  4^^  pair. 


Figure  12.  Minimum  RCRLB  of  the  target  veloeity  [dBm/see]. 


Therefore,  the  performances  of  each  bistatic  system  are 
strongly  related  to  the  configuration  of  the  bistatic  triangle, 
that  is,  to  the  positions  of  the  transmitter,  the  receiver  and  the 
target.  Using  different  transmitting  and  receiving  systems,  the 
target  can  be  seen  by  different  bistatic  configurations. 


Then,  knowing  the  coordinates  of  each  transmitter  and 
each  receiver  of  the  whole  system,  it  is  possible  to  find,  for 
each  point  of  the  analyzed  area,  the  TX-RX  pair  having  the 
best  performance,  in  terms  of  the  lowest  RCRLB. 

Figures  10  and  13  show,  in  a  colour  coded  map,  the 
transmitter-receiver  pair  which  has  the  minimum  RCRLB  for 
each  point  of  the  analyzed  area.  The  color-map  of  these 
figures  is  divided  into  4  colours,  each  of  which  is  associated 
with  one  of  the  4  analyzed  bistatic  systems.  Figures  9  and  1 1 
show  the  minimum  RCRLB  of  the  target  range  and  the  target 
velocity,  respectively,  that  is  the  value  of  the  RCRLB  which 
is  provided  by  the  transmitter-receiver  pair  which  has  the 
minimum  RCRLB. 


V.  Conclusions 

In  this  paper,  we  exploited  the  relation  between  the  AF  and 
the  CRLB  to  calculate  the  bistatic  CRLBs  of  target  range  and 
velocity.  The  bistatic  CRLBs  provide  a  local  measure  of  the 
estimation  accuracy  of  these  parameters.  Then,  we  used  the 
information  gained  through  the  calculation  of  the  bistatic 
CRLBs  for  the  choice  of  the  optimum  transmit-receive  pair  in 
a  multistatic  radar  system.  The  performance  of  each  bistatic 
channel  heavily  depends  upon  the  geometry  of  the  scenario 
and  the  position  of  the  target  with  respect  to  each  receiver  and 
transmitter.  We  approach  the  problem  of  optimally  selecting 
the  TX-RX  pair  based  upon  the  CRLB  for  the  bistatic 
geometry  of  each  pair.  The  optimal  pair  was  defined  as  that 
exhibiting  the  lowest  bistatic  CRLB  for  the  target  velocity  or 
range.  These  results  can  be  used  for  the  dynamical  selection  of 
the  TX-RX  signals  for  the  tracking  of  a  radar  target  moving 
along  a  trajectory  in  a  multistatic  scenario.  Ongoing  research 
focuses  on  dynamic  reconfiguration  of  selected  TX-RX  as  a 
function  of  RCRLB  s  as  well  as  tracking  accuracy.  The  results 
reported  here  have  been  obtained  considering  a  scenario 
composed  of  Passive  Coherent  Radars  which  exploit  a  single 
non  FM  commercial  radio  station  as  transmitter  of 
opportunity. 
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